f 


I 

I 




CSDL-T-926 

PREDICTIVE MOMENTUM MANAGEMENT 
FOR A SPACE STATION MEASUREMENT 
AND COMPUTATION REQUIREMENTS 

by 

John Carl Adams 
August 1986 



JN!lSfA-CR-172026) PREDICTIVE HOHENTOM N88-10866 

HARKGEMERT FOP A SPACE STATION HEASOBEMSNT 
AND COHPUTATIOH 8EQOIHEHBNTS (Draper 

(Charles Stark) Lab.) 148 p Avail; NTIS Onclas 

HC A07/HF A01 CSCL 22B 33/18 0104960 




Predictive Momentum Management for a Space Station 
Measurement and Computation Requirements 


by 

John Carl Adams X 

B.S., Massachusetts Institute of Technology, 1984 


Submitted in Partial Fulfillment 
of the Requirements for the 
Degree of 


MASTER OF SCIENCE 
in 

AERONAUTICS AND ASTRONAUTICS 
at the 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
August 1986 

(^John Carl Adams, 1986 


The author hereby grants to M.I.T. and the C.S. Draper Laboratory, Inc. 
permission to reproduce and to distribute copies of this thesis document 
in whole or in part. 


Signature of Author. 


Approved by. 


Certified by. 


Accepted by. 


(/ Department of 




Aeronautics and Astronautics 
August, 1986 




Michael A. Paluszek 
Technical Supervisor, CSDL 


Prof. Walter M. Hollister 
Thesis Supervisor 



Prof. Harold Y. Wachman 
Chairman, Departmental Graduate Committee 



The Charles Stark Draper Laboratory, Inc. 
Cambridge, Massachusetts 021^9 



2 



ABSTRACT 




Predictive Momentum Management for a Space Station 
Measurement and Computation Requirements 


Submitted to the Massachusetts Institute of Technology 
in Partial Fulfillment of the Requirements of 
Master's of Science 
August 1986 
by John Carl Adams 


An analysis is made of the effects of errors and uncertainties in 
the predicting of disturbance torques on the peak momentum buildup on a 
space station. 

Models of the disturbance torques acting on a space station in low 
earth orbit are presented, to estimate how accurately they can be pre- 
dicted. An analysis of the torque and momentum buildup about the pitch 
axis of the Dual Keel space station configuration is formulated, and a 
derivation of the Average Torque Equilibrium Attitude (ATEA) is pre- 
sented, for the case of no MRMS (Mobile Remote Manipulation System) 
motion, Y vehicle axis MRMS motion, and Z vehicle axis MRMS motion. 

Results showed the peak momentum buildup to be approximately 20000 
N-m-s and to be relatively insensitive to errors in the predicting tor- 
que models, for Z axis motion of the MRMS. The peak disturbance momen- 
tum for no motion and Y axis motion of the MRMS was found to vary 
significantly with model errors, but not exceed a value of approximately 
15000 N-m-s for the Y axis MRMS motion with 1 deg attitude hold error. 

Minimum peak disturbance momentum was found not to occur at the ATEA 
angle, but at a slightly smaller angle. However, this minimum peak 
momentum attitude was found to produce significant disturbance momentum 
at the end of the predicting time interval. 

Thesis Supervisor: Prof. Walter M. Hollister 

Professor of Aeronautical and Astronaut! cal Engineering 

Technical Supervisor: Michael A. Paluszek 

Technical Staff, Charles Stark Draper Laboratory 
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1.0 INTRODUCTION 


NASA's plans for a permanently orbiting space station have raised 
many interesting questions in the area of attitude control, and partic- 
ularly in the area of momentum management. Over a long operational 
lifetime, even the relatively small torques of the environment in low 
earth orbit can provide significant momentum buildup that must be dealt 
with in some manner; either through the use of magnetic torquers, ChG's 
(Control Moment Gyros) , or thrusters. 

There are many sources of unwanted momentum on the space station. 
The two most significant result from the aerodynami'c torques and gravity 
gradient torques acting on the spacecraft. But there are many other 
sources, such as, crew motion, docking, and HRMS (Nobile Remote Manipu- 
lation System) movement. 

The usual method for dealing with this unwanted momentum would be to 
allow the momentum to increase to a certain point and then dump it using 
the RCS jets, magnetic torquers, or just tilting the spacecraft and 
using the gravity gradient torques to counteract this momentum. It has 
been shown [5] though, that the disturbance torques expected on the 
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space station will require a large momentum storage if momentum exchange 
devices are to be used to neutralize these disturbances. 

Because of this, it has been proposed that a new system for control- 
ling the buildup of unwanted momentum be implemented, namely one which 
predicts in advance the disturbance torques on the space station, rather 
than dealing with their effects after-the-fact. Such a predictive 
momentum management system would allow a more optimal placement of the 
attitude of the space station so that disturbance torques might cancel 
each other. This, in general, reduces the momentum storage requirements 
placed upon the momentum exchange devices. 

One example of such a predictive momentum management scheme is to 
fly the space station at an "average torque equilibrium attitude". It 
is desired that the space stations attitude remain constant, but for a 
given period of time what should this attitude be? The average torque 
equilibrium attitude is such that if we can predict the disturbance tor- 
ques on the space station for a time T in the future, then at the end of 
that time the integral of all the disturbance torques, and thus the net 
momentum, will be zero. Flying the space station at such an attitude 
has been shown to reduce the momentum storage requirements for attitude 
control by a factor of 4 (see fig. 1). 



The goal is to find an attitude in which to orient the spacecraft 
such that, for a given time, the peak momentum storage requirement on 
the attitude control system is minimized. 

The problem with a predictive system is that accurate knowledge of 
the torques expected on the spacecraft must be had in advance in order 
to choose an attitude such that the unwanted momentum is zero after a 
given period of time. 

Here is a list of possible sources of disturbance torques; 

1. Gravity Gradient forces 

2. Aerodynamic drag 

3. Docking 

4. HRHS motion 

5. Crew motion 

6. Venting 

7. Solar radiation pressure 

8. Radiation pressure from re-radiated and scattered 

radiation from earth 

9. Frictional torques between rotating and non-rotating components 

10. Changes in inertia due fuel consumption, solar 

panel motion, crew motion, NRMS motion, etc. 

The prediction of these torques requires accurate modeling or meas- 
urement of many of the characteristics of the space station and its 
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environment. Some of tKe more important quantities that need to be spe- 
cified are; 

1. Vehicle Inertia Hatrix- An accurate model which includes the 
changes due to mass shifts involved with HRHS motion, solar panel 
rotation, docking, crew motion, venting and fuel consumption. 


2. Vehicle Drag Coefficient- A model which includes changes due to 
orientation, solar panel motion, shadowing, changes in atmospheric 
condi tions (i .e. composition, density), and configuration 
changes (i.e. HRHS position, docked vehicles) 


3. Vehicle Center of Hass and Center of Pressure- Again, a model 
which accounts for changes in configuration due to docking, HRHS 
activity, etc. 

4. Vehicle Environment- Atmospheric density and composition. Gravita- 
tional anamolies. Earth's radiation and re-radiation, solar radi- 
ation and its fluctuations. 

5. Vehicle Position, Velocity, and Attitude 


16 


2.0 DISTURB AMCE TORQUE MODELS 


This chapter gives an assessment of the sources of disturbance 
momentum on the space station. It presents for each disturbance torque 
the mechanics of the disturbing phenomena, the mathematical models 
involved with prediction of these torques, and the assumptions and sim- 
plifications that have gone into each model. It also provides the 
intended sources for each of the parameters in the models, be it direct 
measurement or estimation, and an application of these torque predicting 
models to an example configuration of the space station; the dual keel 
conf i guration. 
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Figure 1. Space StatioRt Dual Keel configuration 


There are many different disturbance torques that act on the space 
station, some of which can be characterized as forces acting at a cer- 
tain moment arm from the center of mass of the station (r X F) , and oth- 
ers which are more easily character i zed as a change in the space 
station's angular momentum (dH/dt) . In table 1, the disturbance torques 
that will be dealt with in this paper are separated into these two cate- 
gories and into two sub categories. Those torques characterized as 
forces acting at moment arms are separated into those that act as body 
forces and those that act as surface, contact forces. And of the tor- 
ques that are characterized as a change in momentum, some are the result 
of changes in the space station's inertia, and some are due to changes 
in the angular velocity of the station's orbit. 
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dH/dt 


Body 

Forces 


Contact 

Forces 


dl/dt 


dw/dt 


Gravi ty 
Gradient 


Magnetic 


Aerodynamic 

Torques 

Solar 

Radiation 

Pressure 

Earth 

Radiation 

Pressure 

Venting 

Torques 

Docking 

Torques 

Friction 

Torques 

Crew motion 


Docking Non- 

circular 

MRMS motion orbits 

Solar 

Panel 

Notion 

Crew Motion 

Radiator 

Motion 

Consumable 

Depletion 

Antennae 
Motion 

Station 
Growth 


Station 
Reconf i g- 
uration 


The torque, and corresponding angular momentum, disturbance associ- 
ated with each of these phenomena are to be expressed in a coordinate 
frame fixed in the vehicle. Altogether, there are five coordinate 
frames that are of importance to the dynamics of the space station. 

[ 

Their definition and the transformation matrices between them allow for 

I 

the characterization and resolution of the forces and torques that act 
upon the space station. These frames are; 
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Inertial 


Sun-f i xed 

Local Vertical/Local Horizontal 

Vehicle C.O.H. 

Component C.O.H. 

The inertial frame is important because the radiators are to be 
fixed inertiaiiy, so to determine the inertia change due to their 
motion, a transformation is needed from inertial coordinates to vehicle 
coordinates . 

The same reasoning applies to the solar panels, which are to remain 
sun-fixed. A transformation between vehicle frame and a sun-fixed frame 
is needed to specify their motion. 

The basic approach of this anaiysis is that the space station is 
assumed to be a combination of several components, each of which can be 
characterized by its own mass, inertias, and body-fixed coordinate 
frame. In general, a torque or force can be defined in the components 
frame more easily than in the vehicle frame. So in the final analysis 
these component torques must be combined and transformed into the total 
vehicle frame. 
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2.1 GRAVITY GRADIENT TORQUE 


One of the most significant torques on the space station will come 
from the "gravity gradient" forces, which arise due to the fact that 
only the center of mass of an object in orbit is in force equilibrium 
(see Figure 2 on page 21). The net force on any other incremental mass 
is; 



Fq 6 " Z4iiRo±WLrim 


( 1 ) 



( 2 ) 


Mge - y; R X 

|R +RP 


* Radius from the center of the earth to space station 
center of mass (C.O.M.) 

■ Radius from space station C.O.N. to incremental mass dm 

■ earth's gravitational constant (GM^) 

■ Gravity Gradient moment about C.O.M. of station 


(R„2+2(R^.R)+r2)3/2 


(from Law of Cosines) 


R„3 (1 +2 (R^ • R) /R„= + R2/R^2) 3/2 


1/Ro®[1-3(R„-R)/Ro= + .... ] 


(from Binomial Theorem) 


Hqg • 4 R X C-k/Ro®( 1-3(R„*R)/R3 (R„+R)]dm 


0 '"o' ' o 


and since R X (R^+R) ■ R X R^ 


Hgg » Z4 i.[ 1-3(R^-R)/^2] (R X R^)dm 

R..^ 


y; -n/R^s (R X R^)dm 

( -*> 0, since ^Rdm ■ 0 about the C.O.H.) 


+ 4 3^/R^5 (R-R^) (R X R^)dm 
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( 10 ) 


Hqq - 3n2/R„2 4 (R-R^) (R X R„)dm 
where; n ■ 4i7R^ 

y Perpendicular 



Figure 3. Transformation from body frame to LVLH frame 

In order to obtain the gravity gradient torque in the body coordi- 
nates it is necessary to transform the earth radius vector into body 
coordinates. Assuming roil(^), pitch(e)« and yawC*!*) as the three euler 
angles, the transformation matrix from the LVLH frame to the body center 


of mass frame is; 



' LVLH-b 


cos'l'cose s i nes i n4>cos4' + sin'l'cos4> 
|-sin'l<cose -s i nes i n4>s i n't! + cos'l»cos<t> 
si no “Sin<l>cose 


•s i n0cos4'cos<l' + sin<i>sin'l' 
s i nOcos<i>s i n't* + sin4>cos4' 
cosOcos<i> 


( 11 ) 

and using the small angles approximation: cosO » i, » 0 


T 


LVLH-t> 


■ 1 

e4> + '1' 

-e + 4><l< ” 

-<l/ 

-e<M' + 1 

e* + <l> 

e 

-4> 

1 


(12) 


'•’LVLH-b * 


1 ♦ -e 

-<|i 1 

e -4> 


;] 


(13) 


'OLVLH 


0 

0 

R, 


" ^'•’LVLH-b^ •^OLVLH 


Rob " Ro 


-e 

1 


(14) 


X 

y 

z 


R • R ■ (-ex + 4>y + z) R, 


( 15 ) 

(16) 


R X R„ - R„ 


y - <l >2 
-X - ez 
4>x + ey 


(17) 
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(R-R^) (R X RJ « 

-0xy+<ty^+y2+®4>X2-<l>^y2-4>z^ 1 R^^ 

(18) 

ex^“<t>xy-x2+©^x2-e4>y2-e2^ I 


-64>x^+4»2xyH>x2-e^xy+e<l>y^+©y2 J 


flS 

-exy + yz + <l>(y^- 2 ^) *| R^^ dm 



-<l>xy - X2 + e (x^"2^) 1 
«t>x 2 + eyz J 

(19) 

and with the definitions for moments and products of inertia; 


■ ^(y^+22)dm 



ly • y^(x^+z2)dm 

I„ - /^(X2)dm 

(20) 

Ij ■ 7j,(x2+y2)dm 

lyz ■ ^(y2)dm 



Sc the linearized gravity gradient torque in body coordinates becomes; 


Hqq ■ 3n^ 


-*I.y - 


( 21 ) 
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2.2 AERODYNAHIC TORQUES 


The density of the atmosphere falls off exponentially with altitude, 
but even at the space station's nominal altitude of 250 nautical miles 
(450 km) the effects of atmospheric drag can be felt, especially for 
large spacecraft over extended periods of time. However, in this region 
of the atmosphere, the exosphere, the density is so low that the princi- 
ples of contimuum aerodynamics no longer apply. Since the mean free 
path of the molecules at this altitude is greater than the character- 
istic length of the vehicle, each particle interaction with the surface 
of the spacecraft must be considered as independent and uninfluenced by 
another molecule's interaction. Another assumption that can be made is 
that the velocity of the vehicle is much greater than the thermal veloc- 
ity of the molecules. Because of this, the atmosphere can be modeled as 
a molecular beam, with all incoming particle velocities parallel. 
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Mean Atmospheric Temperature as a Function of Altitude 

Figure 4. Structure of Upper Atmosphere 

When a particle collides with a surface, the resulting interaction 
can range from a totally inelastic collision where the particle is 
absorbed by the surface, to a totally elastic collision where the normal 
momentum of the particle with respect to the surface is reversed and its 
tangential momentum is left unchanged, (see Figure 5} 



Figure 5. Specular and Diffuse reflection 


To characterize this range of interactions for a given surface of a 
given material it is customary to define two 'momentum accommodation 
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coefficients', one for the tangential momentum and one for the normal 


momentum. 


o * Tangential momentum « 

accommodation coefficient (22) 

o' ■ Normal momentum ■ p^^ - p^^^ 

accommodation coefficient Pn" ” (23) 


where: 

Pj^ ■ Tangential momentum incoming particles 
Pjp ■ Tangential momentum reflected particles 
p^^ ■ Normal momentum incoming particles 
p^^ * Normal momentum reflected particles 
“ Normal momentum re-emitted particles 

When both of these coefficients are equal to 1, then the particle 
interactions are completely inelastic, which is diffuse reflection. 
When both of the coefficients are 0, then the particle interactions are 
completely elastic, which is specular reflection. In reality they will 
be somewhere in between, with specular reflection generally increasing 
for an increasing angle of incidence between incoming particles and the 
normal to the surface. Some experimental values for o and o' are shown 
i n F igure 6. [1] 
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Figure 6. Some Experimentally Determined Values for o and o' 

The collision process is a complex physical and chemical inter- 
action, and the behavior of absorbed and reflected particles varies for 
each specific surface composition and the type of incoming particles. 
However, a simplifying assumption can be made due to the fact that in 
the steady state, the space craft surfaces will not be 'clean', but will 
be contaminated by a layer of atmospheric particles both stuck and chem- 
ically bonded to the surface. The most notable of these contaminants is 
monatomic oxygen. This layer of contaminants produces a more homogene- 
ous surface condition on the overall spacecraft, and makes the defi- 
nition of the momentum accommodation coefficients easier. 

In general, a layer of surface contaminants greatly raises the tan- 
gential momentum accommodation coefficient, to about the .85 to .90 
region, for all types of surface materials, except at very high angles 
of incidence ('grazing' angles). The normal momentum accommodation 
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coefficient is less affected by surface contamination but this coeffi 
cient tends to vary less from material to material. 


Once the net momentum transfer from particle to surface is charac- 
terized by these momentum accommodation coefficients, then a drag coef- 
ficient equivalent to that for continuum flows can be formulated as a 
function of these coefficients and the surface's attitude with respect 
to the incoming flow. [2] 



x,y,z fixed to surface 
X normal and into surface 


Figure 7. Frame of reference for surface 
The drag coefficients for each direction in the surface are; 
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- C((2-c')Su/^ + (c'/2)v/VT)e*<S‘'-)A2 

+ (2-c ' ) (1 /2+ (Su^) 2) + ( (c/2) /T^ir/T) (1 +er f (Su J ) ] /$= 

Cpy •- (cru/S2) [S/;;?e*<S“x >*+ u^S2 (1+erf (Su^) ) ] 

Cpj - (<tu,/$2) [S/ 7 e*<Su, )* + u^s2 (1+erf (Su^) ) ] (24) 

where u ■ 

■ unit velocity vector in surface frame 
S “ ratio of vehicle speed to average molecular speed («10) 
T^/T ■ ratio of wall temp, to ambiant temp. (*^.25) 



Once these drag coefficients have been calculated for a surface, 
then the forces exerted on that surface follow as; 


^ - 1/2 pV^ACp 
- 1/2 pV^AC 
j - 1/2 pV^AC 


where, 


p * atmospheric density 
A * area of surface normal to flow 
V “ magnitude of velocity 


(25) 


And the torque about the center of mass due to the aerodynamic force 
on each surface is simply; 


■A1 


r, X F, 


(26) 
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To find the total aerodynamic moment it is necessary to sum the 
moments from all the individual surfaces. What complicates this is the 
fact that for different orientations of the vehicle, some parts of the 
spacecraft may be partially or totally shielded from the flow by other 
parts of the spacecraft. This is known as 'shadowing'. Realizing this, 
the total aerodynamic moment about the center of mass is; 

- 2 ; r^ X (27) 

where i • space station surfaces exposed to the flow 

It may be more convenient, however, to define the moment in terms of 
two other quantities, r^p, and C^, where. 


I r * radius vector from center of mass to center of pressure 

Cp > total space station drag coeff. 

- I, 

I 
I 

where, again, i • spacecraft surfaces not shadowed 
Using these definitions, equation (27) can be rearranged to form, 

"a “ '■pp A 1/2 pV^A^Cp (30) 


(28) 

(29) 
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There are three limiting factors in the degree of accuracy with 
which the aerodynamic moment can be predicted. The first limitation 
is in the ability to determine Cp, the total drag coefficient of the 
space station. This is accomplished computationally by simplifying com- 
plex spacecraft geometries into more basic surfaces such as flat plates, 
spheres, rectangular solids, cylinders, etc., and then predicting which 
surfaces will be shadowed and to what extent by using the geometry of 
the space station's configuration and attitude. 

Also, the values for the momentum accommodation coefficients are 
derived experimentally, and the conditions of the space station's 
external surfaces in the space environment will not be an exact match 
for those in an experimental environment. Simplifications in the space 
station’s geometry and the experimental nature of the accommodation coef- 
ficients both limit the accuracy of the total vehicle coefficient of 
drag. 

One way of improving tha accuracy of the momentum accommodation is 
to perform the measurements similar to those done in laboratory molecu- 
lar beam experiments on board the space station itself [34], using the 
incoming atmospheric particles as the molecular beam. This would avoid 
the problem of simulating the surface contamination the spacecraft would 
encounter in orbit, in the laboratory. 
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Another limiting factor is the determination of the spacecraft's 
velocity. This is essentially a navigation problem and it doesn't 
appear that there will be any problem achieving the desired accuracy. 
Using the Global Positioning System, at the proposed altitude of the 
space station, accuracies 45 to 60 feet in position and .2 feet/sec in 
velocity have been shown to be feasible. [37] 

The strongest limitation on the ability to determine the aerodynamic 
torque is in the prediction of the atmospheric density. At the space 
station's altitude the density is extremely variable. In general, den- 
sity decreases exponentially with altitude, but it has also been shown 
to be dependent on the incoming flux of solar radiation, (see Figure 8) 
particularly at the 10.7 cm wavelength, and is also dependent on the 
variations in the earth's magnetic field, (see Figure 8) It can change 
by up to three orders of magnitude over the entire solar cycle. [13] 
Some of the variations are random in nature while others appear regular- 
ly. Some of the regular variations are; 

1. An 11 year solar activity cycle variation 

2. An annual variation 

3. A semi-annual variation 

4. A 27 day solar rotation variation 

5. A diurnal variation due to earth's rotation 

6. Hagnetic disturbance variations 
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ATMOSPHERIC DENSITY (gm/cm^) 

DAYTIME MAXIMUM 

NIGHTTIME MINIMUM 


Figure 9. Density vs. Altitude for maximum and minimum solar 
activity and diurnal variation 

The annual variation in atmospheric density is due to the changing 
composition of the constituent gasses of the upper atmosphere with lati- 
tude and season. 
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The semi-annual variation is a product of the interaction of the 


earth's magnetic field with the solar wind. This effect reaches its 
minimum in July, its maximum in October or November, and a secondary 
minimum and maximum in January and April, respectively. 

The variation with the highest frequency is the diurnal variation. 
The diurnal density 'bulge' reaches its maximum at a point lagging the 
subsolar point by about 30”, or 2 hours. It varies roughly sinusoidally 
I with longitude, about an average value; 

i 

P - + -5 cos(x - 30*)) (31) 

I 

I 

where x ■ the angular separation between the space station and the sub- 

I 

I solar point. 


2.3 SOLAR RADIATION TORQUE 


The incoming solar radiation carries with it some incoming momentum. 
If the radiation is thought of as a stream of incoming particles, pho- 
tons, each with a momentum proportional to its energy, then the situ- 
ation is analogous to the aerodynamic drag case, where incoming 
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particles interact with the projected area of of the surfaces of the 
spacecraft resulting in a net exchange of momentum. 

The analogy is carried further when the types of possible inter- 
actions are considered. As in the case of aerodynamic particle inter- 
actions, photons can be either absorbed, reflected diffusely, or 
reflected specularly, (see Figure 10 on page 37) 



AasonmoN specular reflection diffuse reflection 

Abson>(ion and Reflection of Incident Radiation 

Figure 10. Possible photon interactions with space station 
surfaces 

To characterize what portion of the incident radiation experiences 
each of the different interactions, three coefficients are defined; 

Oj ■ coefficient of absorption 

opp ■ coefficient of diffuse reflection 

<T|jj • coefficient of specular reflection 


where 

® ‘'a»®ro'®rs 

< 1 

and 

®RD ®RS 

- 1 


These coefficients represent the fraction of the incoming momentum 
that is ei ther ’ absorbed, reflected diffusely, or reflected specularly. 
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I 


I 

I 

for a given surface of a given material, and like the momentum accommo* 
dation coefficients, they are determined experimentally. [26] 

The total energy flux radiated from the sun is known as the 'solar 
constant', S, and is » 1.35X10^ J/m^-sec at the radius of the earth's 
orbit. The solar constant varies annually by about 6 % due to the eccen- 
tricity of earth's orbit around the sun. The momentum flux is; 

i P • S/c (32) 


p is in the outward radial direction from the sun. and c is the speed 

of light, « 3X10® m/sec. 

i 

I 

If we define a coordinate frame that is fixed to the surface of the 
spacecraft; (see Figure 11 on page 38) 

incoming 

, momentum 


Figure 11. Surface coordinate frame, incoming momentum 



I 

I 
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with being a unit normal into the surface, and with Up being a unit 
vector in the incident momentum direction, then the force imparted to 
the i*^ surface by each type of interaction is; 

specular reflection- ■ 2pcose (Acose) u^ 

• 2S/c cos^eAu^ (33) 

diffuse reflection- ■ p(Acose) (Up + (2/3) u^) 

• S/c (Acose) (Up + 2/3 u^) (34) 

absorption- F^^ ■ p (Acose) Up 

» S/c (Acose) Up (35) 

Specular reflection is simply a reversal of the normal component of 
the incoming momentum, while diffuse reflection causes the net momentum 
transfer to be along the the direction (Up + 2/3 u^) . Absorption causes 
a force along the direction of incoming momentum, Up 

The total force on this surface can now be expressed by incorporat- 
ing the absorption and reflection coefficients: 
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■ Solar radiation pressure force 


* * ‘'fiS^RSI ‘^RD^RDi 

■ SAcose/c C(oa+ord)Up + 2 (op5Cose+o„p/3) u J (36) 


Expressing the force in a coordinate frame that is normal and tan- 
gential to the surface can be accomplished by replacing the vector Up by 
Up ■ coseu^ * sineu^. The resulting expression is; 

■ SAcose/c [ (a^+Ogp) coseu^ 

+ (cr^sine+agp(2/3 +sine)+2op5cose)u^] (37) 


The associated moment for that surface about the center of mass is. 


H 


R1 


n X Fg, 


(38) 


And summing across the entire vehicle gives the total moment due to 
radiation pressure: 


^R " ^ * ^Ri 


(39) 
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where the summation is over those surfaces exposed to the sunlight, thus 
eliminating 'shadowing' effects. 


2.4 EARTH EMITTED AND SCATTERED RADIATION TORQUES 


The sun is not the only source of radiation which can place a torque 
on the space station. Another important source of radiation is the 
earth itself. The earth can radiate in two ways, first as a black body 
it emits in the infrared range, and second it scatters incoming solar 
radiation. On average, the earth scatters 34t of the incoming solar 
flux[9], and if we assume that the earth eventually re-radiates all of 
the energy that is incident upon it, the amount of radiation energy 
emitted in the IR range is 66% of the incoming solar energy flux. The 
corresponding momentum flux from this radiation is; 


Psc 


(40) 


PiR • p^R^!jL6§) 


(41) 
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I 


where; 


Rg ■ Radius of the earth 

Rg ■ Radius from center of earth to space station 
■ Incoming solar momentum flux 


The incoming solar radiation is intercepted by the projected area 
the earth presents to the sun, vR,^. The next assumption is that both 
the scattered and emitted radiation are radial in direction with the IR 
radiation being emitted over the entire surface of the planet, and the 
scattered radiation only being emitted from the illuminated side of the 
planet. The amount of momentum reaching the spacecraft must therefore 
be divided by the surface area into which the radiation is being 
divided, which is 4irR^^ and 2irR^^ respectively. The assumption that the 
scattered momentum is radial becomes suspect near the division between 
day and night in the orbit. 

The forces and torques due to th i s' momentum flux from the earth can 
be found by using the same procedure as outlined for the solar radiation 
torque. The momentum flux defined above is simply substituted for that 
of the direct solar radiation. But already it can be seen that these 
torques will be much smaller than those from direct solar radiation, due 
to the division by the surface area of a sphere with the orbit's radius. 
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2.5 MAGNETIC TORQUES 


When in orbit, there will be an interaction between the earth's mag- 
netic field and the net magnetic dipole of the space station, arising 
from current loops or permanent magnetism. Although the spacecraft may 
be designed such that the net magnetic dipole is approximately zero, 
there will be some residual that will interact with the earth's field to 
produce a torque: 


N^ - D X B 


(42) 


where ■ magnetic torque about c.o.m. 

0 “ residual magnetic dipole of space station 
B “ earth's magnetic field (in LEO «.1 gauss) 


To accurately determine this torque it is necessary to be able to 
accurately define both D and B The earth's magnetic field can be mod- 
eled as a magnetic dipole centered in the earth, with a dipole strength 
with an axis that protrudes from the earth's surface at 78.5* N lat- 
itude, 69.0* W longitude. It varies with altitude as 1/R^ and is also 
time varying due to the bombardment of charged particles from the solar 
wind. The field is also latitude and longitude dependent, but still a 
complete model can be formulated. 
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The magnetic field potential is; 


V(r,e,<l>) ■ (a/r)"*^ 2^,,.^" (g^"’cosm4> + h„"sinm<}>)P„"(o) (43) 

where; 

a ■ equatorial radius of the earth 

“ Gaussian coefficients (empirical) 
r ■ geocentric distance 
e ■ colatitude 
4> • east longi tude 
Pf,"'(e) ■ Legendre functions 

The n>1 terms are called 'dipole', the n-2 'quadrupole' , and so on. 

It is easier and more accurate in view of the time variations due to 
"magnetic storms", to directly monitor the magnetic field vector with an 
on-board magnetometer. 

Determining the magnetic dipole of the space station is not as sim- 
ple a measurement. Dipoles arise from three sources; 

1. Current loops and permanent magnetism 

2. Eddy currents 

3. Hysteresis effects 

The first source is the dominant one, and also the one most easily 
accounted for in design, either in an attempt to null the dipole or use 
it to produce a control torque. What is needed is a value for the resi- 


dual magnetic dipole for each component of the space station, so that a 
value for the total residual magnetic dipole can be obtained by adding 
the component values vector iaily. Whether these dipoles can be measured 
in orbit or must be measured at the time of the components manufacture 
is unclear, as is the question whether they will vary significantly with 
time. 

The average expected magnetic torque on the dual keel space station 
is on the order of .0001 N-m, while the maximum expected torque, assum- 
ing a failure which causes the main y axis truss to be the current loop, 
is on the order of .01 N-m. 


2.6 INERTIA CHANGE TORQUES 


The desired attitude for the space station is one which keeps the 
stations z axis aligned with local vertical. This will require that the 
station maintain a constant angular velocity equal to the spacecraft's 
orbital angular veloci ty (w^) . The angular momentum of the space craft 
will be: 

H - Iwg (44) 
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Any change in the space stations mass or configuration win lead to 
a change in the total inertia matrix of the vehicle. Since the stations 
angular velocity is constrained to remain at orbit rate» the angular 
momentum must change, and thus a disturbance torque results on the 
spacecraft. 

+ (dI'7dt)o.p + I"do.p/dt 
dt 

There are several ways in which the space station's inertia can 
change. These are some, ranked according to how big an impact they have 
on the total inertia; 

1. Station growth - module addition 

2. Station reorganization - module movement 

3. Docking - Shuttle orb iter or OTV 

4. Motion of the Remote Manipulation System (RMS) 

5. Solar panel motion 

6. Radiator motion 

7. Crew motion 

8. Consumable depletion - fuel .water, air.food.garbage.etc. 

9. Antenna, and other platform, motion 

10. Motion due to structural flexibility 


(45) 

(46) 
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When the space station is considered as a combination of rigid body 
components, each of which can be characterized by its own inertia and 
mass, then any change in the total vehicle inertia can be described as 
some combination of these four sources; 


1. Component Rotation 

2. Component Translation 

3. Change in Component Inertia 

4. Change in Component Mass 


The inertia of the entire space station can be found as a function 
of the component inertias, the component masses, m^, and the radi- 
us vector from the overall vehicle center of mass to the component cen- 
ter of mass, Rg. 



Figure 12. Relationship of C.O.H. Component and C.O.M. Station 
Frames 


I ' ■ T I" T^ 

A c 'c-b^ c' c-t 


( 47 ) 
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where I''c * component inertia matrix in component frame 
I'jj ■ component inertia rotated into vehicle frame 
and ■ rotation matrix from component frame to vehicle frame. 



“ separation between x^ and x^ axes 
‘ ^ R'c, ♦ K'c 

1, - ^ ♦ R*.. 

(49) 

X. - ^ ♦ R“e,' 

The component moments of inertia translated into the vehicle axes 

are; 




- 

*C2 


^ * CX ^ 

T iCy 

^ CZ ^ Z 


(from Parallel Axis Theorem) 


(50) 


And the component products of inertia translated into vehicle axes 

are; 


I ■ I' 

I 

*cxy * 

I - I’ 

,cxy 

•^cxz ^ 

cxz 

I “I' 

I 

•^cyz * 

cyz 

Ic - Ic' 

+ II 


+ m^Rj^Ry 
+ m^Rj^Rj 
+ m^RyRj 


m^ 


[ 


»v 

•Rc^Rgy 

*^cx^cz 


R2 


cz 


*^cx^cy 

*^cy^cz 


cz 


“^cx^cz 


cy J 


(51) 


(52) 


which can be written more compactly as> 
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le - I'c + in,((R,-R^)£ - R^RJ) 
where £ is the identity matrix 


( 53 ) 


Total station inertia: 


I - StdV'TT + m^((R^-R^)£ - R,RJ)) 


( 54 ) 


In order to find the inertia change torques, this expression for the 
vehicle's inertia matrix must be differentiated with respect to time. 


dl/dt - \ [d/dt(Tl^TT)+ dmg/dt((R^*R^)£ - R^Rj) 

+ m^d/dt ( (R^ • R^) i - RcR J) ] (55) 

- 2; Ed/dt(T)ITT + Td/dt(Ig)TT + TIgd/dt(TT) 

+ dmg/dt((R^-Rg)£ - R^RJ) 

+ m^(d/dt(Rg-Rg)£ - d/dt(R^Rj)) (56) 

For simplicity, ' has been shown simply as I^, and the transfor- 
mation matrix, T^.^, has been shown simply as T. 

And since d/dt(Rg*Rg) - (57) 

where ■ velocity of component center of mass with respect to the cen- 
ter of mass of the vehicle. 

And d/dt (R^RJ) - V^Rj + R^Vj (58) 

The derivative of the inertia matrix in the vehicle frame becomes; 
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(59) 


dl/dt - [d/dt(T)ITT + Td/dtac)T'^ + Tl^d/dt (T^) 

■ + dn»c/dt{(Rc-Rc)i - 

+ m^{(2V^-R^)£ - V,RJ - R^Vj) 


So there are eight parameters which need to be specified to deter- 
mine rate of change of the inertia of the space station so that the tor- 
que due to these changes can be found. The/ are; 


Vb* «l(T^.^)/dt 
Rc* V, 

Ic. 

fflg, dm^/dt 


(Component Rotation) 

(Component Translation) 

(Change in component inertia) 
(Change in component mass) 


2i6.1 Solar Panel Torques 

As an example of the torques that can be produced on the space sta- 
tion due to a time varying moment of inertia, here is an analysis of the 
most predominant motion, the once per orbit rotation of the solar pan- 
els. The figure shows the coordinate frames and relative motion with 
respect to the rest of the space station for this motion. 
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Figure 13. Relative motion of solar panels 


As a component of the entire space station, the solar panels are 
constrained such that their center of mass does not translate with 
respect to the spacecraft. Their mass and inertia are constant, and so 
equation (59) reduces to 


dl'Vdt - (dC/dt)JC^ + CJ(dC/dt)T 


(60) 


The torque on the vehicle due to this changing inertia is; 


M - VI' '“o + (‘*1' '/dt)«o (61) 

- VI"“o + C(dC/dt)JCT + CJ(dC/dt)T]«^ (62) 


dC/dt "F-sinede/dt 0 cosede/dt" 
0 0 0 
_-cosede/dt 0 -sinede/d^ 


(63) 


where: 
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e “ ftt 
do/dt * ft 


dC/dt ■ -ft Psine 

0 

-cose” 

0 

0 

0 

cose 

0 

sine_ 


(dC/dt)^ ■ -ftPsine 0 cose” 
0 0 0 
^cose 0 sine 


The solar panel inertia matrix, J. is diagonal 


( 64 ) 


(65) 


J - Fj, 0 0 “ 

0 0 
0 0 


( 66 ) 


So the torque becomes; 


H « u)/i' + -ft r . . . 0 • • • 1 r ° ”1 

0 I I — n I 

I ... 0 ...Jo 


( 67 ) 


With the vehicle angular rate only about the y axis, the second term 
in the torque equation is zero and only the euler coupling term is left. 


0 - 
0 
0 


"1 

oj L^13- 


I ** 
‘*‘13 


I " 

■*■12 
T II T II 

■*■22 ‘*‘23 

T II T II 

‘*‘23 ■*‘33 


0 

-n 

0 


( 68 ) 


53 



( 69 ) 


H 


n^I " 

n XJ3 
0 

-n^I " 

n 


The torque due to the solar panel rotation is the same as if they 
weren't rotating, that is simply the euler coupling torque from the 
vehicles rotation at orbit rate. This is because the cross terms I. 12 " 
and Ijg" in the total vehicle inertia matrix are not dependent on the 
motion of the solar panels. It is also interesting to note that there 
is no pitch torque, about the y axis, due to the motion. If friction is 
included, there will be a pitch torque due to the solar panel motion. 


2.6.2 Toroue due to HUMS motion 

Determining the inertia change due to HRHS motion requires a few 
assumptions. First, the HRHS itself is assumed to be a point mass with 
no rotational inertia. The mass of the HRMS is assumed to be constant, 
and its motion is constrained to the y-2 plane in space station coordi- 
nates, along paths that are parallel to either the y or 2 axes. 

The equation for the change in inertia of the vehicle (59) reduces 
to 


dl/dt - m^p^s((2v„*trji - v„rj - r^vj) 


(70) 
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where and equal the HRHS position and velocity with respect to 
the core/solar panels center of mass. The possible MRHS velocities con- 
sidering the constraints on the motion are; 


V 


M 



(for y-motion) 



(for 2 -motion) 


And the HRMS position is of the form; 


r 


n 



(71) 


(72) 


(73) 


The rate of change of the inertia matrix for motion of the HRHS par- 
allel to the y and 2 axes becomes, respectively; 


dl/dt 


"'mrms 


2r 

0 

0 


V 

my^my 


-r 


0 

0 


• r w 

?? "V 


(74) 


dl/dt 




2r 


V 

nz^iiz 

0 

0 


0 


2r 

-r 


V 

nz^mz 

V 

ny'^mz 



(75) 
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And the torques due to this HRHS motion are, 


^MRMS “ 


^MRMS " 


where; 



0 

2"*»»RMs'’mz''m2'’ 

”'"«Ras'”ay''a2^ 


( 76 ) 


(77) 



There is no pitch torque produced 


(78) 

(79) 

(80) 

by the NRHS motion along the y 


axis, but there is a torque about the 2 axis. The motion along the z 


axis produces torques about both the y and 2 axes. 


A plot of the pitch torque produced for an HRHS manuever when the 
space station is held stationary in LVLH is shown in Figure 14 The 
maneuver is motion along the 2 axis of the spacecraft from z«+20m to 
z*-20m with a speed of .5 m/s. 
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Figure 14. HRHS motion torques, z axis maneuver 


2.7 FRICTIONAL TORQUES 


Since the space station core is to remain nominally earth pointing, 
while the solar panels are sun pointing, there will have to be some kind 
of rotating hinge at the joint between the two. This interface will 
have to transmit the power that is being generated by the solar panels 
to the core and will necessarily involve some friction between the two. 
This friction will cause an internal torque on the space station that if 
correctly modeled could be useful in the decreasing unwanted disturbance 
momentum. 
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There are two types of power transmission that are likely to be used 
at solar panel/core interface. One is slip rings which transmit DC 
power and the other is rotary transformers which transmit AC power. In 
the case of slip rings, friction arises due to the carbon copper contact 
across which the current passes, while in the case of rotary transfor* 
mers, the friction occurs in the bearings which support the rotating 
coils of the transformer. 

In either case, a simple model of the frictional forces is that they 
are proportional to the normal force applied between the two contacting 
surfaces. The constant of proportionality is the coefficient of kinetic 
friction, which is approximately constant for varying values of 
solar panel angular velocity. 

Fp - ( 81 ) 
The torque associated with this frictional force is; 

Hp - r X Fp (82) 

where r is the moment arm at which the friction force is applied. It is 
necessary to know the geometry of the contact at the interface in order 
to determine this moment arm. It is also necessary to know the geometry 
of the contact in order to determine in what direction normal forces are 
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being applied. The net force on the hinge will be the net acceleration 
of the space station at that point, times some mass, as yet to be deter* 
m i ned . 

In order to determine the friction torque, it is necessary to know 
the coefficient of friction, the geometry of the contact, and the accel* 
eration of the spacecraft at the hingepoint. The structure and size of 
the slip rings or rotary transformers are constrained by the amount of 
power they are meant to transfer. The total power production on the 
space station is on the order of 150 kw, which is divided between two 
hinges. The interfaces are sized for roughly 75 kw of power, and this 
wilt determine the diameter of the slip rings or transformers and thus 
the moment arm of the friction torque. The total torque will be the sum 
of the torques for each hinge. 

A simplifying assumption can be made that the acceleration of the 
space station at the hingepoint will be small, so the normal forces will 
be small and can be assumed to be constant. The friction in these types 
of devices, as determined from an ESA study by the European Space Tri- 
bology Laboratory [36] is on the order of 1 to 10 N-m. 
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2.8 SUMMARY 


Here are the maximuni predicted levels of disturbance torque for all 
the phenomena described in this paper; 

Table 2. Maximum Disturbance Torques 


TORQUE 

MAXIMUM VALUE (N-m) 

Gravity Gradient 

21 .93 

Aerodynamic 

.119 

Solar Radiation 

.023 

Earth Radiation 

.052 

Magnetic 

.01 

MRMS Z-Motion 

202.47 

Solar Panel Friction 

10.0 
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a.O ANALYSI S OF HQHEMTUH BUILDUP; PITCH TORQUES ALONE 


In order to predict the momentum response to disturbance torques of 
the entire vehicle, it is useful to first model the single axis response 
of the space station. Since the pitch axis wilt experience the largest 
aerodynamic torques, due to the offset of the center of pressure created 
by the OTV and OMV berths in the dual keel space station configuration, 
this is the axis that will be modeled. 

Once modeled, the momentum response of the vehicle to these torques 
can be found by integrating them over time, and the peak momentum the 
disturbance torques will produce can be found from examining this func- 
tion for momentum. This peak momentum will be the determining factor in 
the sizing of the momentum exchange devices that will be needed to deal 
with this disturbance momentum. Once the peak momentum is determined, 
it is then possible to predict the effect that various attitudes have 
upon the peak momentum value, the effect that uncertainties in attitude 
have upon the peak value, and the^effect that uncertainties in the tor- 
que models have on the peak momentum value. 
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3.1 SIMPLE AERO/GRAVITY GRADIENT MODEL 


The disturbance torques being considering in this section are; 

1. Aerodynamic 

2. Gravity Gradient 

3. MRMS Motion 


The total disturbance torque then becomes; 

H[3 “ M^ + Mgg + Nmfus (83) 

Mg - Mg(e.t) (84) 

The pitch torques can be modeled with differing levels of complexi- 
ty. depending on whether effects such as the change in inertia of the 
space station with solar panel and MRMS motion are to be included. A set 
of simple torque models follows; 

M^ » -1/2 PoV^Cg(A,pr,p + A^r^) cosed + 1/2 cos(nt)) (85) 

Hqg “ 3n2[(i^ - IJsinecose - 1^^] (86) 

"mrms " 2M^r„^v„^n - 2M^nv„^ (r„„ + v^^t) (87) 

(for 2-motion of MRMS) 


where: 


62 



e ■ pi tch angle 

t * elapsed time in orbit from point of maximum atmospheric density 
n ■ orbit rate ■ 2w/T 

B average atmospheric density 

V • magnitude of spacecraft velocity with respect to atmosphere 

Agp ■ projected area of the solar panels 

Ag * projected area of the core 

Cp ■ overall vehicle drag coefficient 

r^ B core aerodynamic moment arm 

r,p B solar panel aerodynamic moment arm 

I B vehicle inertia 

"*MRMS * >nass of MRUS 

r„, B position of MRMS with respect to C.O.M. 
v„ B velocity of MRMS with respect to C.O.H. 


The assumptions made in formulating these simple equations for the 
disturbance torques on the space station are; 

1. Solar Panels not rotating 

2. Distance to C.O.P. constant over time 

3. The drag coefficient doesn't change with attitude 

4. The inertia of the vehicle doesn't change with solar panel and 
HRMS motion 

5. There are no shadowing effects in the aerodynamic or radiation 
torques 

6. The atmosperic density variation is a simple cosine model of the 
diurnal component of the variation 


The total momentum built up over time is simply the integral of 
these torques over time. 


H(e,t) - Ho(e,t) dt 


(88) 
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If we consider only the aerodynamic and gravity gradient torques. 


and then only the simple models for these torques, the expression for 
the momentum of the spacecraft as a function of time and attitude can be 
found. 

(89) 

(90) 

Where; 

(91) 

(92) 

H(e,t) - dt 

■ .^^[K,cose(1 +1/2 cos(nt)) + K,sinecose - K-Jdt 

(93) 

H(e,t) ■ (K^cose)t + (K^2n cose) sin (nt) + (Kjsinecose) t - Kgt (94) 

H(e,t) • (K^cose + K^sinecose - Kj)t + (K,/2n cose) sin (nt) (95) 

For this simple model, the disturbance momentum is just the superpo- 
sition of a linear and a sinusoidal term. 


K, - -1/2 PoV^Cp (A,pr,p + A^r^) 
Kj - 3n2(I. - I^) 

K, . 3n^„ 


■ K,cose(1 +1/2 cos(nt)) 
Mg- ■ K-sinecose - K- 
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The Average Torque Equilibrium Attitude (ATEA) is, for pitch alone, 
the angle at which if the spacecraft is fixed for a given time T, the 
net momentum buildup is aero. 

H( 04 ,T) ■ (K,cose^)T+ (K,/2n cose^)sin(nT) + (K 2 Sine^cosej)T - KgT « 0 

(96) 

H(e^,T) ■ (K^cose^ + KgSine^cose^ - Kg)T + (K^2n eoseJsin(nT) (97) 

If we choose the time T to be a multiple of the orbit period T, 
where T*2w/n, then the sin(nT) term is aero, and we are left with this 
condition to solve for e^; 


K^cose^ + Kjslne^cose^ - Kg ■ 0 (98) 

This is a transcendental equation, and as such cannot be solved 
explicitly for the ATEA, e^. But if the constants K^, Kg, and Kg are 
known, then an iterative solution for can be set up. 

e* - sin'’(1/Kg(K3/cose^ - K,)) (99) 

To find the peak momentum over the time period T, it is necessary to 
examine the extremums of the function H(e,t) by taking the derivative of 


H(e,t) with respect to time and setting it equal to zero. Differentiat- 
ing with respect to time yields the original torque expression. 

dH/dt ■ [K^cose (1 +1/2 cos(ntp)) + K^sinOcose - Kg] 

“ (K^cose + KgSinecose - Kg) + (K^/2 cose) cos (ntj ■ 0 

( 100 ) 

t ■ 1/n cos"* (K g - K ^ cose - K^sinecose) ■ time of peak momentum 

(K,/2)cose (101) 

By substituting this expression for the time of peak momentum build- 
up on the space station into the expression for the momentum (95), the 
peak momentum buildup as a function of the attitude at which the space- 
craft is held can be found. 

Hp “ (K^cose + KgSinecose - Kg) (1 /n) cos'* (K g - K ^ cose - Kgsinecose) 

K^/2 cose 

+ 1/n / ( (K,/2) cose) ^ - (Kg - K^cose - Kgsinecose)^ (102) 

It is now possible to examine the effect of keeping the space sta- 
tion held at various attitudes on the value of the peak momentum accumu- 
lated over one orbit. The first attitude to be examined is the LVLH 
hold attitude. In this attitude, e * 0, and the equation for peak 
momentum reduces to; 
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Hp * (K, - Kj) (1/n)cos-’(2(K3 - K^)/K,) + (1/2n)^K,2 - 4(K3 - 

(103) 

This value exists only if jKg - Kj ^ K^/2 which means that the 
cosine terms in the momentum equation are not overpowered by the linear 
terms in such a way that there are no extremums of the function in the 
intervai 0 S t ^ T. Even if the vaiue of exists, it may stiii oniy 
be a iocai maximum or minimum, and intuitiveiy. the maximum momentum 
bui idup for one orbit shouid be at the end of the orbit, at t ■ T. At 
the end the momentum is; 

Hp - (K, - K 3 )T + (K,/2n)sin(nT) 

- (K, - Kg)T (104) 

If the space station is heid in ATEA. then the peak momentum 
expression reduces to; 

Hp • K,/2n cose^ (105) 

To compare these two peak vaiues. it is necessary to evaluate the 
constants numerically. The parameters for the duai keei space station 
are defined in Tabie 3 on page 68. 
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Table 3. Dual Keel Parameters, Solar Panels Fixed, No MRMS 


Ajp - 1974.64 m2 

r^p ■ +.3709 m 

Ag - 261 .17 m2 

r^ ■ "5.0672 m 

I„ - 62105091 kg-m2 

I, - 45995243 kg-m2 

I„ - 5930594 kg-m2 

■ 10000.0 kg (w/payload) 

- 167869.2017 kg 

N,p - 5510.1999 kg 

Pp ■ 1.5 X 10'''2 kg/m® 

V ■ 7624 m/s 

n - 1 .1188 X 10'® s*'' 


For these values the ATEA pitch angle is; 
e* - -23.6* 

And the peak momenta are; 

Hp(LVLH) - -124585.64 N-m-s 
Hp(ATEA) - 35.28 N-m-s 

For this simple model, the ATEA hold is much better that the LVLH 
hold. This is because the Dual Keel is an essentially aerodynamical 1y 
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balanced design, so the dominant torques come from the gravity gradient. 
The ATEA is such that the gravity gradient torques are in equilibrium, 
so when the space station is moved away from ATEA and held in LVLH, a 
significant disturbance momentum is created that must be dealt with. 
That the difference is so large is only the case of this simple model, 
which is subject to the assumptions outlined earlier. 


a. 1.1 Simple Model with Solar Panel Rotation 

If the solar panel motion is considered, the torques are of the same 
form as before, except that now the projected area of the solar panels 
to the incoming atmospheric particles is a function of time. 


Ma - ■ 

-1/2 {A,pr,p + Aprp) cosed + 1/2 cos(nt)) 

(106) 

Mqg " 

3n^[(Ij - I^)sin0cose - 1,^^] 

(107) 

where 




A,p |cos(nt)| 

(108) 


So the torque becomes; 
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Mj, * K^cose (1 +1/2 cos (nt) ) 

+ K/ cosed +1/2 cos(nt)) / cos (nt) , for 0 £ t i v/2n 

<-cos(nt) , for ir/2n S t 5 3ir/2(109) 

( cos (nt) , for 3w/2n i t £ 2it/n 

+ K2sinecose - Kg 


Where; 



- -1/2 pv^CpA^rg 

(110) 


' - -1/2 pv2CpA,pr,p 

(111) 

Kg 

■ 3n2(Ij - Ij^) 

(112) 

^3 

- 3n2l^, 

(113) 


Integrating this expression for torque over time yields a new 
expression for the disturbance momentum that takes into account the 
rotation of the solar panels. 


H(e, t) ■ K^cose (t + 1/2n sin(nt)) 

+ K,'cose/(1/n sin(nt) + t/4 + 1/8n sin(2nt)) 

\ 0 £ t ^ w/2n 

<(2/n + ir/4n - 1/n sin(nt) - t/4 - 1/8n sin(2nt)) 

I */2n S t S 3ir/2n 

( (4/n - ir/2n + 1/n sin (nt) + t/4 + 1 /8n s i n (2nt) ) 

3w/2n S t i 2ir/n 

+ (Kgsinecose - Kg) t 

(114) 
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The assumption that the variation in spacecraft inertia due to the 
rotation of the solar panels can be ignored is valid for the pitch axis. 
The variation of the vehicle's inertia about this axis is zero since the 
rotation of the solar panels is constrained to be about this axis. How- 
ever, the variation of the spacecraft's inertia about the other two axes 
due to solar panel motion can be significant, and since the pitch gravi- 
ty gradient torque is dependent on these off axis inertias, the assump- 
tion must be examined more closely. The total vehicle inertia is; 


I - I' + CJCT - Kr/r/ (115) 

The time varying portion of this inertia is the CJC''’ term, which 
expands to, 


CJC^ 


J^cos^(ftt) + JgSin^(nt) 0 

0 Jj 

(Jg - J,)sin(ftt) 0 


(Jg - J ,) s i n (nt) cos (fit) 

0 (116) 
J,sin^(flt) + Jgcos^ (fit) 


The gravity gradient torque is a function of the I^, 1^, and 1^^ 
terms of the inertia matrix. 


- I,,' + J,cos2(nt) + JgSin2(fit) - I^' + J, + (Jg 
Ij ■ Ij' + J,cos2(nt) + JgSin2(nt) ■ Ij' + Ji + (Jg 
^**3 " Ji)sin(nt)cos(m) 


J^)sin2(nt(117) 
J^) cos^ (nt (118) 
(119) 


The inertias are all dependent on the difference between the Jg and 
terms of the solar panel inertia. The closer these two values are, 
the less impact the solar panel motion has on the gravity gradient tor- 
que. 


Some predicted values for solar panel inertia are; 

- 16545372.76 N-m^ 

Jg - 1292047.280 N-m^ 

1/ - 45559718.78 N-m^ 

Ij' - 30679373.76 N-m® 

- 5930594.86 N-m^ 

Jj - J, - -.25 

J, + I,' (120) 

■— “ ”*32 

J, + I,' (121) 

and 

Jg - J, - -2.57.32 

I,.' (122) 

The ratios show that the time varying components of the inertias are 
significant in comparison to the inertia of the vehicle if the solar 
panels are considered as non-rotating. 
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3.2 COHPLEX TORQUE MODELS 


A more complete modeling of the torques on the space station would 
take into account the inertia changes associated with the motion of the 
solar panels and HRMS, the resulting changes in the moment arms of the 
aerodynamic and radiation torques due to these mass shifts, and the tor- 
ques due to the changing inertia of the space station as it is being 
held at a constant rate. 

The motion of the solar panels produces no net torque about the 
vehicle y-axis, but it does cause the 1^, and 1^^ vehicle inertias 
to be time varying, which affects the gravity gradient torque. Similar- 
ly, motion of the HRMS along directions parallel to the pitch axis pro- 
duces no pitch torque, but does contribute to the changes in the gravity 
gradient torques. The motion of the HRMS parallel to the 2 -axis has the 
most pronounced effect on the spacecraft. It produces a sizable torque 
as well as both changing the vehicles inertia and changing the pitch 
moment arm of the aerodynamic torque by shifting the center of mass of 
the entire vehicle parallel to the z-axis. (F igure 15 on page 74) 

The MRMS maneuvers are shown in Figure 16 on page 75. 
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-38.4 


MRMS « 10000 I K« 


Figure 15. Shift of C.O.H. with MRMS motion 
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fO. -17» 88. -38.4) 

Z-Motlon Start 



(0, -17.88,+54.a^ 


♦ Z-Motlon Stop 
4 (0, -17.88, +54.86) 


Y-Motion Stop Y-Motion Start 


Figure 16. NRHS maneuvers 


3.2.1 Compl ex Model Y-Metion of MRMS 


The aerodynamic torque, for y-motion of the NRHS, is the same as for 
the case of no motion of the HRHS, except that the moment arms of the 
pitch torque are adjusted for whatever initial z position the NRHS has, 
which in this case is *•■54.86 m. 


M* ■ -1/2 (A^r^ + A,pr,p) (1 +1/2 cos (nt)) cose 


(123) 



where; 


Agp * frontal area of the solar panels 

> A,p| cos (nt) I , due to solar panel rotation 

The aero torque can be expressed as a function of time and the pitch 
angle. 

■ K^cose (1 +1/2 cos(nt)) 

+ K/ cosed +1/2 cos(nt)) I cos (nt) , for 0 ^ t S ir/2n 

<-cos(nt) , for ir/2n S t S 3ir/2(124) 

( cos (nt) , for 3'w/2n S t S 2it/n 

For simplicity, the parameters in the equation have been collected 
into the constants and K/ 

(125) 

(126) 

The p'itch gravity gradient torque is a function of the I,* and 
vehicle inertias. 

Hge - 3n2[(i^(t) - I^(t))sinecose - I«(t)] (127) 


K, - -1/2 pv^CoAprp 

and 

K,' ■ -1/2 pv^C|jA,pr,p 
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To determine the torque, the inertia must be found as a function of 


time, given the motion of the NRMS and the solar panels. 


I - I' + CJC'f - Kr^’^rT* 


(128) 


where; 


K - 

and 


(129) 


r^ B distance from core/solar panel center of mass to HRHS 
J > solar panel inertia matrix ■ constant 
C * transformation from solar panel to core coordinates 
I' B core inertia matrix, plus the inertia due to the mass 
offset from solar panels ■ constant 

C “ C (t) and r^ “ r^(t) 


Again, the time varying portions of the inertia expression are; 


CJC^ ■ I J^cos^(nt) + JgSin^(nt) 

0 

(Jg - J^) si n (nt) cos (nt) 


(Jg - J^) sin (nt) cos (nt) 
0 

J^sin^ + Jgcos^ (nt) 


and 


(130) 


Kr,^ * - K 


-(r 


Tz^ r,/) 
'-Tx''Ty 
’’tx''Tz 


-(r 


- '^Tx'^Ty- 
2 . 1 . r 3 
Tz ^ '^Jx 

'■Ty*'Tz 


•’Tx’'tz 

" (^Ty 


Tx 


(131) 
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For the y-motion of the MRMSj * 0 

'■ti “ '■tzo “ constant 

•“Ty • ’’Tyo + 


This gives the time varying inertia components; 


Ix(t) ■ I^' + J^cos2(nt) + JgS indent) + K (r^^^ + (-132) 

Ij (t) ■ Ij' + J^sin^(nt) + J,cos^ (nt) + K (r^ (133) 

I„(t) ■ !„' (^3 " J,) sin (nt) cos (nt) (134) 


And the gravity gradient torque is; 


Mqq ■ (Kj- + Kgcos (2nt) ) s inecose ~ Ks “ ^ 3/2 sin(2nt) (135) 


Again, for simplicity, the parameters have been collected into the 
constants Kg, Kg, K^, and Kg, where the constants are defined as fol- 
lows; 


Kg 

K 4 

•<5 


3n2(Jg - 
3n^ (Krij- 
3n2 (I ■) 



(136) 

(137) 

(138) 

(139) 


Adding the aerodynamic torque to the gravity gradient torque gives 
the total torque on the space station, for y-motion of the HRHS, as a 
function of time and pitch angle. Integrating that expression over time 
yields the disturbance momentum imparted to the spacecraft that will 
have to be dealt with by the CHG's. 
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H(o, t) * K^cose (t + 1/2n sin(nt)) 

+ K/cose|(1/n sin(nt) + t/4 + 1/8n sin(2nt)) 

\ 0 i t S ir/2n 

\(2/n + It/An - 1/n sin(nt) - t/4 - 1/8n sin(2nt)) 

I ir/2n $ t S 3‘ir/2n 

\(4/n - •ir/2n + 1/n sin(nt) + t/4 + 1/8n sin(2nt)) 

3ir/2n S t S 2ir/n 

+ ((Kj- K4) t + Kg/2n sin(2nt)) sinecose - Kgt + Kg/4n cos (2nt) - Kg/4n 

(140) 

This expression can be seen to reduce to the simpler expression for 
momentum that didn't take into account inertia changes from the internal 
motions, by setting the Kg and constants to zero. 

Now that we have this expression for the disturbance momentum, the 
questions that need to be answered are; 

1. What is the ATEA (Average Torque Equilibrium Attitude) at which 
to fly the space station to cancel this momentum over a given 
time T? 

2. What is the peak momentum expected for various attitudes? 

3. What effects do uncertainties in the attitude and in the torque 
models have upon this peak momentum prediction? 
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The peak momentum is important as it will determine the sizing of 
the momentum exchange devices that will be used to compensate for this 
disturbance momentum. 

The ATEA pitch angle can be found by employing its definition; 

H(e^, T) • 0 , T » 2‘ir/n ■ orbit period (141) 

K^cose^(T) + K/cose^(4/n) + (K^- K 4 )Tsine^cose^ - KgT - 0 (142) 

Solving this transcendental equation for can be accomplished 
iteratively by using this expression; 

- sin'‘‘[1/(K2~ K 4 ) ((Kg/cose^) - (K^ + K/ 2/ir))] (143) 

The time t*T was chosen because it simplifies the equation for 
however any time could have been used, while multiples of T also give a 
simpler express i on. 


3 . 2.2 Complex Model Z-Wotion of MRMS 

For motion of the MRMS parallel to the 2 -*axis, different expressions 
for the torque and momentum must be obtained. The variations of the 
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moment arms of the aerodynamic torque with a shifting center of mass, 
and the torque due to a changing must be included along with the var- 
iation due to changing 1^, and inertias. 


The aerodynamic torque for z-motion of the MRMS is; 

« (K, - K/t) cosed +1/2 cos(nt)) 

( + (Kj - Kj't) cosed +1/2 cos(nt)) cos (nt) 

for 0 £ t !£ t^ 


( (K, - K^'t,) cosed +1/2 cos(nt)) 

+ (Kj " Kj't,) cosed + 1/2 cos(nt))^ cos (nt) ,for t^ S t S ir/2n 

<-cos (nt) ,for w/2n S t i 3ir/2n 
( cos (nt) ,for 3ir/2n S t S 2ir/n 

for t, i t i 2ir/n (144) 

Where; 

K, - -1/2 PoV^CpA^r^^ (145) 

K/ - -1/2 PeV^CoA^(V^-^%p)''Mz 

Kj - -1/2 PoV2CoA,pr,p„ (147) 

Kj' - -1/2 PoV2CeA,p(M^/Mp+M,p)v^, (148) 

and t^ ■ time the NRHS motion stops 


The gravity gradient torque for MRUS z-motion becomes; 
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Mqg * 1 (^ 3 - - Kgt^ + Kgcos (2nt)) sinecose - K, - Kg/2 sin(2nt) 

I for 0 S t S 

*(Ka- K^t, - Kgt ^2 + Kgcos (2nt) ) sinecose - K, - Kg/2 sin(2nt) 

for t, S t S 2ir/n (149) 

Where; 


Kg - 3n2(I^' - I,' - KrT„2) (150) 

K 4 « 3n2 (2KrTjeVnj) (151) 

Kg - 3n2(Kv„j2) (^52) 

Kg - 3n2{Jg - j^) . (153) 

K, - 3n2(I^^') (154) 

The torque from the changing due to the' HRMS motion is; 

^MRMS * (r^j^Viuj + (155) 

(for 2 -motion of HRMS) 


Integrating the sum of these three torques yields the expression for 
the disturbance momentum for the 2 -motion of the NRHS; 
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H(e, t) ■ K,cose (t + 1/2n sin(nt)) 


I 


+ KjCOsed/n sin(nt) + t/4 + 1/8n sin(2nt}) 

+ K/cose(t^/2 + 1/2n^ cos (nt) + t/2n sin(nt) - 1/2n^) 

+ Kj'coseCt/n sin(nt) + 1/n^ cos (nt) + t/8n sln(2nt) 

+ 1/16n2cos(2nt) + tV8 - 17/16n2) 

+ (Kgt - K4t^/2 - Kjt®/3 + Kg/2n sin(2nt)) sinecose 
- Kgt + Kg/4n cos (2nt) - Kg/4n 

+ 2M^n (r.,^v„t + v„V2 t^) 

for 0 S t i t, 


(K, - K/t,)cose(t + 1/2n sin(nt)) 

^ Q 2 

+ {(Kg - K^t, - K5t,2)t + K-/2n sin (2nt)) sinecose 
- Kgt + Kg/4n cos (2nt) (K4t//2 + 2/3 Kjt,®) sinecose - Kg/4n 

+ 2M„n (r^,^v„t, + v^^/2 t,2) 

+ (Kj - Kg' t^) cose/ (1/n sin(nt) + t/4+ 1/8n sin(2nt)) 

\ for t, S t S ir/2n 

<(2/n + w/4n - 1/n sin(nt) - t/4 - 1/8n sin(2nt)) 
) for ir/2n S t S 3it/2n 

I(4/n - ir/2n + 1/n sin(nt) + t/4 + 1/8n sin(2nt)) 
for 3ir/2n S t S 2w/n 

for t^ £ t £ 2ir/n (156) 

where 

Qg • K,'cose(t,V2 - 1/2n^ cos (ntg) - 1/2n^) (157) 

+ Kg'cose(t,V8 - 1/n^ cos(nt^) - 1/16n^ cos(2nt^) + 17/16n^) 




And lastly, the ATEA angle for this maneuver can be found by an 
iterative solution of an equation of the form; 
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» sin'^l/C^ (Cg/cose^ - (C, + Cj + C 3 ) ) ] (158) 

where; 

C, - (K, - K/t^)T (159) 

Cj - (Kj - K2't^)4/n (160) 

C 3 ■ Qj/cose^ (161) 

C 4 - (K 3 - - Kgt,2)T + (K4t,2/2 + 2/3 (162) 

Cg - -K,T + 2M„n(r,„v„,t, + v„,V2 t,=) (163) 
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4.0 RESULTS FROM MODELS 


In this section, the torques and momentum response from the simple 
model with solar panel rotation are presented in Figure 17 on page 86 to 
Figure 24 on page 93 . 

The torques and momentum response from the complex model with no 
motion of the MRHS are presented in Figure 25 on page 94 to Figure 33 on 
page 102 . 

The torques and momentum response from the complex model with Y*mo~ 
tion of the MRHS are presented in Figure 34 on page 103 to Figure 42 on 
page 111 . 

And the torques and momentum response from the complex model with 
Z-motion of the MRHS are presented in Figure 43 on page 112 to Figure 51 
on page 120. 
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Figure 18. Homentutn, Simple Hodel, LVLH hold 




Figure 19. Momentum, Simple Model, ATEA - 1 deg hold 
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Figure 24. Torque, Simple Hodel, LVLH hold 
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Figure 30. 
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Figure 31. Momentum, No Motion MRMS, ATEA hold, p 
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Figure 32. Torque, No Notion HRHS, LVLH hold 
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Figure 34. Homentum. Y-Hotion MRHS, LVLH hold 
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Figure 37. Momentum, Y-Motion HRMS, ATE A - 1 deg hold 
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Figure 39. Momentum, Y-Motlon MRMS, ATEA hold, p + 501; 
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Figure 40. Momentum, Y-Hotion MRMS, ATEA hold, p - 50)1 
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Figure 41. Torque, Y-Hotlon HRMS, LVLH hold 
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Figure 42. Torque, Y-Hotion NRMS, ATE A hold 
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Figure 43. Nomentum, Z-Notion MRNS» LVLH hold 
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Figure 44. ' Homentinn, Z-Motion HRNS» ATEA hold 
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Figure 45. Momentunit Z-Motlon KRNS, Min. Peak Momentum Attitude 
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Figure 46. Momentum, Z-Motion MRMS, ATE A - 1 deg hold 
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Figure 47. Momentum, Z-Motion HRMS, ATEA + 1 deg hold 
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Figure 48. Homentum, Z>Hotion HRNS, ATE A ho id, p + 50% 
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Figure 49. Homentum, 2-Motion MRMS, ATE A hold, p 
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Figure 51. Torque, Z-Motion MRMS, ATEA hold 


5.0 DISCUSSION OF RESULTS 


5.1 ATTITUDE AND PEAK NON5NTUH 

The peak disturbance momentum varies significantly with the attitude 
at which the space station is held. (Table 4 on page 123) For all the 
models, and all possible motions of the MRMS, It can be seen that it is 
significantly better to hold the spacecraft fixed at the ATEA angle than 
to hold the vehicle axes aligned with the LVLH frame. The factor of 
four reduction observed for the simple model In [5] is now a factor of 
-fifty reduction in the peak disturbance momentum due to. the more aero- 
dynamical ly balanced design of the dual keel space station design when 
compared to the previous power tower design. 

Comparing the peak momentum observed with the simple modei and that 
observed with the complex model with no MRHS motion shows the importance 
of including the effects of the motion of the solar panels on the iner- 
tia of the space station in the analysis. This leads to an increase in 
the expected peak momentum on the order of a factor of fifty. 
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The better performance at the ATEA attitude can also be attributed 
to the fact that the vehicle axes of the dual keel design are not the 
principal axes of the space station. A large xz product of inertia 
exists when the vehicle axes are aligned with the LVLH frame which pro- 
duces a large pitch gravity gradient torque. 

The most important result of examining the effect of attitude on the 
peak disturbance mcxnentum is the discovery that, except for the simple 
torque model, the ATEA attitude does not produce the minimum peak dis- 
turbance momentum on the space station. The ATEA angle is close to the 
attitude which produces the minimum peak momentum (Table 5 on page 124), 
but the actual angle that produces the minimum peak momentum is less 
than the ATEA angle by up to one degree. The reason for this Is that 
while the ATEA attitude guarantees that the disturbance momentum will be 
zero at the end of the time interval T, it does not guarantee that the 
median disturbance momentum will be zero during that interval. The 
angle which makes the median disturbance momentum zero is that which 
produces the minimum peak momentum. . 

Fixing the space station's attitude at the angle which minimizes the 
peak disturbance momentum does not necessarily give zero momentum at the 
end of the time interval. Whatever momentum is left must be dealt with 
on the next time interval, so this minimum peak momentum attitude may 
not be the best attitude. And as can be seen is Table 4 on page 123, 
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the reduction in peak momentum from the ATEA hold case is only 10-20% 
for a1 1 the models. 


Table 4. Peak Homentum vs. Attitude Specification 


CASE 

Peak Homentum 
(N-m-s) 

% Difference 
From Nominal 

Simple model, ATEA hold 

44.1927 


Simple model, LVLH hold 

-124503.0 

**■281627.53% 

No HRNS motion, ATEA hold 

2370.39 


No MRNS motion, Ninimum Peak 

Momentum Attitude 

2010.40 

-15.2% 

No HRHS motion, LVLH hold 

-124503.0 

•*■5152.43% 

Y motion MRHS, ATEA hold 

2467.73 


Y motion MRMS, Minimum Peak 

Momentum Attitude 

2040.70 

-17.3% 

Y motion MRMS, LVLH hold 

-121480.0 

•*■4822.74% 

Z motion MRMS, ATEA hold 

20445.0 


Z motion MRMS, Minimum Peak 

Momentum Attitude 

18430.0 

-10.3% 

Z motion MRMS, LVLH hold 

BHB 

•*■415.87% 


I 
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Table 5. 


Atti tude 

Pitch Angle 

Simple Model ATEA hold 

-23.57820760 deg 

No MRHS motion ATEA 

-26.27417614 deg 

No MRMS motion min. 
peak momentum attitude 

-25.73385 deg 

Y-motion MRMS ATEA 

-7.64112094 deg 

Y-motion MRMS min. 
peak momentum attitude 

-7.51450 deg 

Z-motion MRMS ATEA 

-6.28786112 deg 

Z-motion MRMS min. 
peak momentum attitude 

-5.18090 deg 


5.2 HRHS MANEUVER AND PEAK MOMENTUM 


The motion of the MRHS increases the peak disturbance momentum in 
general. Motion parallel to the y-axis of the space station produces a 
slight increase in the peak momentum, while motion parallel to the 
2-axis and for the full length of the space station increases the peak 
momentum by an order of magnitude. (Table 6 on page 126) 

It also can be seen that, for 2 motion of the MRHS, the peak momen- 
tum is dependent on the speed at which the MRHS moves. Faster speeds 
decrease the peak momentum slightly. Taking advantage of this has its 
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limits in that the HRhS will have some maximum speed at which it can 
move. Also, the analysis of the momentum buildup shown here does not 
take into account the starting and stopping torques of the HRHS, which 
become more significant as its operating speed is increased. 

Also, if the HRHS 2 motion is less than the entire length of the 
vehicle, then the peak disturbance momentum resulting from the motion is 
greatly reduced. If the distance traveled is only 20 meters instead of 
90 meters, the peak momentum can be reduced by a factor of 10, as seen 
in Table 6 on page 126. 

Because the motion of the HRHS parallel to the 2 axis produces such 
a large peak momentum, it might be more efficient to treat it as an iso- 
lated disturbance and deal with it individually, rather than include it 
in any prediction scheme where an optimal attitude is being chosen. 
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Table 6. Peak Momentum vs. MRMS Maneuver 


CASE 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

No MRMS motion, ATEA hold 

2370.39 


Y motion MRMS, ATEA hold 

2467.73 

+4. Ilk (over 
no motion) 

Z motion MRMS, ATEA hold 

20445.0 

+762. 5k (over 
no motion) 




Z motion. Faster, Vm*.5 m/s 

18791.2 

-8.09k 

Z motion. Slower, Vm*.125 m/s 

23664.7 

+1 5.75k 

Z motion. Shorter distance 

2913.45 

-85.75k 

Z motions Shorter, Slower 

2939.27 

-85.62k 

Z motion. Stop at C.O.M. 

-12728.5 

-37.74k 


5.3 ATTITUDE UNCERTAINTY AND PEAK NONENTUM 


Another question that must be answered after the most desirable 
attitude is determined is how closely can this attitude be followed. A 
reasonable estimate of the accuracy of the control of the space station 
is that it will be able to hold an attitude to within -t-1 deg or -1 deg 
error. The resulting effect on the peak momentum of being that far off 
from the ATEA angle is shown in Table 7 on page 128. 
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It can be seen that the attitude uncertainty has a substantial 
effect on the peak momentum in the no HRHS motion caset and a very large 
effect on the peak momentum in the Y motion case, but the Z motion case 
is relatively insensitive to errors in attitude. 

The effect on the Y motion peak momentum of an error in attitude is 
important because with a relatively small error the peak momentum value 
is raised to nearly that of the Z motion case. This is probably due to 
the large aerodynamic moment arm that is produced when the MRHS is sta- 
tioned at the lower keel, as it is in the y axis manuever. 

There is no longer an advantage of a lower peak momentum when deal- 
ing with the y axis maneuver versus the z axis maneuver due to this 
error in the attitude. 
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Table 7. Peak Momentum vs. Uncertainty in Attitude 


CASE 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

No MRMS motion, ATEA hold 
(pitch angle ■ -26.274176 deg) 

2370.39 


No MRMS motion, ATEA * 1 deg 

-3524.19 

+48.681^ 

No MRMS motion, ATEA - 1 deg 

4605.99 

•••94. 3U 

Y motion MRMS, ATEA hold 
(pitch, angle ■ -7.6411209 deg) 

2467.73 


Y motion MRMS, ATEA + 1 deg 

-15586.6 

•••531 .B2% 

Y motion MRMS, ATEA - 1 deg 

15437.4 

•►525.57% 

Z motion MRMS, ATEA hold 
(pitch angle ■ -6.2878611 deg) 

20445.0 


2 motion MRMS, ATEA 1 deg 

18544.1 

-9.3% 

Z motion MRMS, ATEA - 1 deg 

23244.7 

-►13.69% 


5.4 AERODYNAMIC MODEL UNCERTAINTY AND PEAK MOMENTUM 


The real uncertainty in all of these momentum models is the aerodynamic 
torque and momentum contribution. As already noted, the atmospheric 
density and the vehicle drag coefficient are unpredictable in the short 
term due to effects such as solar activity and shadowing of one part of 
the space station by another. The effect of errors In different aero- 
dynamic parameters in the model is shown in Table 8 on page 129, Table 9 
on page 130, and Table 10 on page 131 
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Variations from the predicted values of the aerodynamic parameters 
have little effect on the no motion MRMS and z motion HRHS peak disturb- 
ance momentum values. The effect on the y motion case is greater, with 
a 50% variation in the atmospheric density producing a 25% increase in 
the peak momentum value. In general, the contribution to the peak 
momentum is small so variations in the aerodynamic model have little 
effect. The y motion case is affected the most because for that case 
the moment arms of the core and solar panels are the largest. 

Table 8. Peak Momentum vs. Uncertainty in Aerodynamic Parameters 


CASE 

No HRHS Notion 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

Actual density ■ pred. density 
(Pred. density - 1.5X10E-12) 

2370.39 


Act. den. ■ Pred. den. - 50% 

2309.69 

-2.56% 

Act. den. « Pred. den. 50% 

2461 .39 

* 3 . 64 % 

Act. drag ■ Pred. drag - 20% 
(Pred. drag coefficient ■ 2.7) 

2345.93 

-1 .03% 

Act. drag ■ Pred. drag ■*’ 20% 

2394.84 

+1 .03% 

Act. SParea"Pred. SParea - 25% 
(Pred. S.P. area • 1974.64) 

2393.43 

+ .97% 

Act. SParea"Pred. SParea 25% 

2347.35 

-.97% 

Act. Carea ■ Pred. Carea - 25% 
(Pred. CORE area • 261.17) 

2316.92 

-2.26% 

Act. Carea * Pred. Carea 25% 

2427.72 

-2.42% 
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Table 9. Peak Momentum vs. Uncertainty in Aerodynamic Parameters 


CASE 

Y Motion MRMS 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

Actual density ■ Pred. density 
(Pred. density ■ 1 .5X10E-12) 

2467.73 


Act. den. ■ Pred. den. - 50% 

1889.96 

-23.41% 

Act. den. ■ Pred. den. 50% 

3051 .58 

••■23.66% 

Act. drag ■ Pred. drag - 20% 
(Pred. drag coefficient ■ 2.7) 

2235.77 

-9.40% 

Act. drag “ Pred. drag + 20% 

2700.62 

••'9.44% 

Act. SParea*Pred. SParea - 25% 
(Pred. S.P. area ■ 1974.64) 

2280.85 

-7.57% 

Act. SParea">Pred . SParea * 25% 

2654.83 

•►7.58% 

Act. Carea “ Pred. Carea - 25% 
(Pred. CORE area - 261.17) 

2364.57 

-4.18% 

Act. Carea * Pred. Carea * 25% 

2571 .59 

•►4.21% 
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Table 10. Peak Momentum vs. Uncertainty in Aerodynamic Parameters 


CASE 

Z Motion MRMS 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

Actual density ■ Pred. density 
(Pred. density ■ 1.5X10E-12) 

20445.0 


Act. den. ■ Pred. den. - 501 

20223.7 

-1 .08% 

Act. den. ■ Pred. den. + 50% 

20666.4 

+1 .08% 

Act. drag “ Pred. drag - 20% 
(Pred. drag coefficient ■ 2.7) 

20356.5 

-.43% 

Act. drag ■ Pred. drag + 20% 

20533.6 

•*■.43% 

Act. SParea"Pred. SParea - 25% 
(Pred. S.P. area ■ 1974.64} 

20441 .7 

-.02% 

Act. SParea*Pred. SParea + 25% 

20448.3 

+ .02% 

Act. Carea ■ Pred. Carea - 25% 
(Pred. CORE area - 261 .17) 

20337.6 

-.53% 

Act. Carea ■ Pred. Carea + 25% 

20552.4 

+ .53% 


5.5 HASS PROPERTIES UNCERTAINTIES AND PEAK NOHENTUH 

In Table 11 on page 133, Table 12 on page 134 , and Table 13 on page 
135 the effects of uncertainties in the values of the masses and iner- 
tias of the space station upon the peak value of the disturbance momen- 
tum are shown. In each ease the spacecraft was held in the ATEA 
attitude. 
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It can be seen that the Z motion case is relatively insensitive to 
changes in either the inertia components or the masses of the different 
elements. None of the cases responds with a significant change in peak 
momentum to a change in the core or solar panel mass, but both the no 
motion and y motion case peak momentum values are increased significant- 
ly by an lOt uncertainty in the inertia components of the core and solar 
panels. 

Another interesting effect that can be seen in the graphs of the 
momentum response of the space station to errors in the inertia values 
of the models, is that the ATE A angle is hot affected by changes in the 
inertia values of the solar panels. This is because their cyclic motion 
produces a net zero momentum and so doesn't require any compensation for 
in the ATEA. 
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Table 11. Peak Momentum vs. Uncertainty in Mass and Inertia 


CASE 

No MRMS Motion 

Peak Momentum 
(N-m-s) 

) Difference 
From Nominal 

Actual Mass SP ■ Pred. Mass SP 
(Pred. Mass SP ■ 5510.1999 kg) 

2370.39 


Act.MassSP “ Pred.MassSP - lOt 

2365.51 

-.21) 

Act.MassSP ■ Pred.MassSp 10); 

2375.24 

4.20) 

Act. MassC • Pred. MassC - 10); 
(Pred. Mass CoreH67869.2 kg) 

2375.78 

+ .23) 

Act. MassC “ Pred. MassC 10) 

2365.95 

-.19) 

Act. I2P - Pred. IZP - 10) 
(Pred. IZP - 30679373. 76kg-m2) 

25681 .3 

+983.42) 

Act. IZP - Pred. IZP + 5) 

-12840.7 

+441 .71) 

Act. IZP - Pred. IZP + 10) 

-25578.3 

+979.08) 

Act. IXP - Pred. IXP - 5) 
(Pred. IXP - 4555971 8. 78kg-m2) 

-19068.7 

+804.45) 

Act. IXP - Pred. IXP + 10) 

38137.4 

+1508.91) 

Act. IXZP • Pred. IXZP - 10) 
(Pred. IXZP - 5930594. 86kg-m2) 

12506.9 

+427.63) 

Act. IXZP - Pred. IXZP + 10) 

12507.0 

+427.63) 

Act. J1 ■ Pred. J1 - 10) 

(Pred. J1 • 16545372.76 kg-m2) 

-855.16 

-63.92) 

Act. J1 - Pred. J1 + 10) 

5531 .05 

+133.34) 

Act. J3 - Pred. J3 - 10) 

(Pred. J1 • 15315869.66 kg-m2) 

5296.18 

+123.43) 

Act. J3 - Pred. J3 + 10) 

-620.27 

-73.83) 
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Table 12. Peak Momentum vs. Uncertainty in Hass and Inertia 


CASE 

Y Motion MRMS 

Peak Momentum 
(N-m-s) 

% Difference 
From Nominal 

Actual Mass M ■ Pred. Mass M 
(Pred. Hass MRMS - 10000 kg) 

2467.73 


Act. Mass M ” Pred. Hass M - 10% 

-7904.86 

••>220.33% 

Act. Mass M ■ Pred. Mass M + 10% 

7819.11 

+216.85% 

Act.MassSP “ Pred.MassSP - 10% 
(Pred. Mass SP ■ 5510.1999 kg) 

2473.18 

+ .22% 

Act.MassSP * Pred.MassSP 10% 

2462.31 

CM 

CM 

• 

1 

Act. HassC ■ Pred. HassC - 10% 
(Pred. Mass Cor e«1 67869. 2 kg) 

2460.52 

-.29% 

Act. HassC * Pred. HassC 10% 

2365.95 

+ .27% 

Act. IIP » Pred. IZP - 10% 
(Pred. IZP - 30679373. 76kg-m2) 

8526.52 

+245.52% 

Act. IZP - Pred. IZP + 10% 

-8526.53 

+245.52% 

Act. IXP - Pred. IXP - 10% 
(Pred. IXP - 4555971 8. 78kg-m2) 

-12662.1 

+413.11% 

Act. IXP » Pred. IXP + 10% 

12662.1 

+413.11% 

Act. IXZP - Pred. IXZP - 10% 
(Pred. IXZP - 5930594. 86kg-m2) 

12506.9 

+406.82% 

Act. IXZP - Pred. IXZP + 10% 

12507.0 

+406.82% 

Act. J1 ■ Pred. J1 - 10% 

(Pred. J1 - 16545372.76 kg-m2) 

-1042.87 

-57.74% 

Act. J1 - Pred. J1 + 10% 

5288.25 

+114.30% 

Act. J3 ■ Pred. J3 - 10% 

(Pred. J1 - 15315869.66 kg-m2) 

5078.50 

+105.80% 

Act. J3 “ Pred. J3 + 10% 

-844.10 

-65.79% 
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Table 13. Peak Momentum vs. Uncertainty in Hass and Inertia 


CASE Peak Momentum % Difference 

Z Motion MRMS (N-m-s) From Nominal 

Actual Mass M ■ Pred. Mass M 20445.0 

(Pred. Mass MRMS - 10000 kg) 

Act. Mass M - Pred. Mass M - 10^ 18047.0 -11. 73^^ 

Act. Mass H - Pred. Mass M + lO^t 22845.2 ,7A% 

Act.MassSP - Pred.MassSP - lOi^ 20468.5 

(Pred. Mass SP 5510.1999 kg) 

Act.MassSP - Pred.MassSP 4 lOt 20421.7 -.lit 

Act. MassC “ Pred. MassC - lOt 20211.5 -1.14t 

(Pred. Mass Core*1 67869.2 kg) 

Act. MassC ■ Pred. MassC + lOt 20640.5 *«-.96t 

Act. IZP - Pred. IZP - 10% 20988.7 +2.66% 

(Pred. IZP - 30679373. 76kg-m2) 

Act. IZP - Pred. IZP + 10% 19966.1 -2.34% 

Act. IXP - Pred. IXP - 10% 19733.7 -3.48% 

(Pred. IXP - 4555971 8. 78kg-m2) 

Act. IXP - Pred. IXP + 10% 21496.4 +5.14% 

Act. IXZP - Pred. IXZP - 10% 21844.1 +6.84% 

(Pred. IXZP - 5930594. 86kg-m2) 

Act. IXZP • Pred. IXZP + 10% 19594.5 -4.16% 

Act. J1 - Pred. J1 - 10% 19740.2 -3.45% 

(Pred. J1 - 16545372.76 kg-m2) 

Act. J1 - Pred. J1 + 10% 22104.4 +8.12% 

Act. J3 - Pred. J3 - 10% 21924.6 +7.24% 

(Pred. J1 - 15315869.66 kg-m2) 

Act. J3 - Pred. J3 + 10% 19792.6 -3.19% 
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5.6 CONCLUSION 


In alt these cases, the Z motion of the NRNS has produced the larg- 
est values of peak momentum for ATEA hold. If the effect of the NRNS 
motion is to be included in the predictive momentum management scheme, 
then the values of peak momentum for the Z motion case will be the val- 
ues used to size the momentum exchange devices necessary to compensate 
for the disturbances. It can be seen that the values of peak momentum 
for the Z motion of the NRNS are relatively insensitive to errors in the 
models or measurements that went into predicting what attitude is best 
for decreasing the peak momentum value. Because of this insensitivity 
it appears that the predictive scheme will work even with the inaccura- 
cies in the models and uncertainties in the measurements. But this is 
at the cost of accepting the highest expected momentum values, and 
being required to carry momentum exchange devices sized accordingly. If 
NRNS motion were compensated for independent of the predictive momentum 
management system, then the peak momentum value could be reduced by a 
factor of five, to that of the expected peak disturbance momentum for 
the no NRNS motion case with a 1 deg attitude uncertainty. 

Also, the results show that the goal of minimizing peak momentum is 
not necessarily the best goal since it leaves a large disturbance momen- 
tum at the end of the predicting time interval. And while the ATEA 
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angle doesn't produce the minimum peak disturbance momentum, it does 
come relatively close, as well as leaving zero disturbance momentum at 
the end time. 


5.7 RECONHENOATIONS FOR FUTURE WORK 

The minimum peak momentum attitude must be defined explicitly and a 
compromise between the ATEA attitude and the minimum peak momentum atti- 
tude be found. 

Also the effect of expanding the time scale of the prediction of 
torques should be examined. How does it affect the uncertainty in the 
environment and how does it affect the peak momentum and ATEA calcu- 
> 1 at ion? 
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APPENDIX A. EQUATIONS OF MOTION 


The model of the space station from which the models in the momentum 
analysis were derived is one which consists of three interconnected 
rigid bodies. These bodies are the core of the space station (habita- 
bility, logistics and laboratory modules as well as the supporting truss 
structure), the rotating solar panels, and the mobile remote manipu- 
lation system (NRHS) . As weli as being rigid, each component in the 
composite structure is assumed to have a constant inertia and mass, in 
its own frame. The NRHS is assumed to be a point mass so its inertia is 
negl igible. 

This model was chosen because it will allow the for the examination 
of the effects of two important internal disturbance torques on the 
attitude motion of the space station; the friction torque between the 
rotating solar panels and the core, and the inertia change torques due 
to the motion of the NRHS, as well as the effects of external torques 
such as the gravity gradient and aerodynamic torques. 

To start the derivation of the equations of motion of the three-body 
space station, first define the angular momentum of the composite space- 
craft, which is the momentum of each component about its own center of 
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mass, plus the moment of the linear momentum of each component about the 
composite center of mass. Bodies 1, 2, and 3, are the core, solar pan- 
els, and NRHS, respectively; 

“ A[I^u) + Ij (w'Kl) + l3«o 

+ m, (r, X dr^dt) + mj (rj X drj/dt) + mg (rg X drg/dt) (164) 


where r.,, rg, and rg are the vector positions of the three bodies with 
respect to the vehicle center of mass. In an inertial frame. 

Since the solar panel center of mass remains fixed with respect to 
the center of mass of the core, these two bodies can be combined into 
one, with the core/solar panels now being body 1 and the NRHS now being 
body 2; 


■ A[(I + CJC''^)ii) + CJn + m^r^** (w*r., + dr,/dt) 

+ mgrg^Cw^rg + drg/dt)] 


(165) 




Figure 52. The definition of 1 and "two" bodies 


where; 

I > inertia of. the core structure about its owm c.o.m., 

plus the inertia due to the mass offset of the core and solar 
panels from their combined center of mass 
J B inertia of the solar arrays about their own c.o.m and fn-their 
own frame 

(i> B angular velocity of the entire space station, in the body frame 
n B angular velocity of the solar arrays with respect to the core 
expressed in the solar panel frame. 

C * transformation matrix from solar panel frame to the core frame 
A B transformation matrix from core frame to the inertial frame 
m., B mass of core and solar panels combined 
B mass of the HRHS 

r., B vector from total space station c.o.m. to the c.o.m. of the 
core and solar panels combined 
r^ B vector from total space station c.o.m. to the HRHS 
(and r^, r^ are now in the body frame) 
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The notation ( ) ** refers to the skew symmetric matrix that is con- 
structed from a column matrix. When multiplied by another vector The 
skew symmetric matrix gives the matrix equivalent of a vector cross pro- 
duct. 

Since the centers of mass of the NRNS, core/solar panels, and the 
total vehicle remain col inear, the whole system can be represented by 
one vector from the core/solar panel center of mass to the NRNS. 



BODY 1, CORE /SOLAR 
PANELS 


Figure 53. The definition of r.^ in terms of r^ and rj 


(m^+mj) 


rj • "IiIt 


(166) 
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Using these definitions for r, and Tj, the momentum equation 
becomes, 

Hg - [(I + CJC^ - KrT*rT*)w + CJn + Kr^’^dr^/dt] (167) 

where 

■ vector from the c.o.m. of the combined core and tolar panels 
to the HRHS 

K ■ m^ffl2/(m^ * mj) ('reduced mass' of system) 

Differentiating the angular momentum with respect to time in the 
core frame yields an expression for the torque on the vehicle in the 
core frame. 

M - I"do)/dt + + «*cjn + C(n*j - Jft*)C^« + CJdft/dt (168) 

+ K (ci)’'r^’‘drT/dt + ry*«*drj/dt 
+ dr^/dt^'w^r^ + r^^d^r^/dt^) 


where; 

I" • I + CJC^ - Kr^^ry* (inertia of total vehicle) 

In differentiating this equations, the following matrix identities 
were used; 
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dA/dt - u)*A 


(169) 


dC/dt - n^c 


(170) 


In order to solve this equations for the angular velocity of the 
space station, expressions are needed for the time rate of change of the 
solar panel rate, dn/dt, and the time rate of change of the HRHS veloci- 
ty, d^Tj/dt^. These expressions are obtained by writing the equations 
of motion of the solar panels and the NRHS separately. The equations of 
motion for the solar panels alone are; 


- ACJ (ft + C'^w) (171) 

Mp - dHp/dt (172) 

Mp • w^CJn + w^CJC^w + CJdft/dt + CJC^dw/dt + (ftxj - Jft*)C^o) (173) 

dft/dt - J-’cT [Mp - C(ft’<J- JfV‘)CV W^CJC^W- CJC^dw/dt] (174) 

where Mp is the disturbance torque plus the frictional torque exerted on 
the solar panels at the hinge between the solar panels and the core of 
the space station. The equations of motion for the NRMS alone are; 


*T “ 


(175) 
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ds^/dt ■ + Adr^/dt 


(176) 


d^Sj/dt^ • A<i)*<D*‘r^+ A(d«/dt)’‘r^ + 2Ao)’‘dry/dt + Ad^r^/dt^ ■ (177) 

where ■ HRHS inertial acceleration 
and Sf * NRMS inertial position vector 

Fn • dVj/dt + + (d«i>/dt) ’^r^ + w^w^^r^ (178) 

where F^ is the acceleration of the NRNS in the body frame. 

d^rj/dt^ « dv.f/dt ■ F„ - 2«*dr^/dt- (dw/dt)*r.f - «*‘«*‘r^ (179) 

The three derived equations of motion, for the KRNS, solar panels 
and for the total vehicle, all assume that there are no constraints on 
the relative motion between the components of the space station. In 
fact, the motions of the HRHS and the solar panels are very constrained. 
The soiar panels can rotate only about the Y axis of the vehicle, and so 
must have aero position, rate and acceleration about the other two axes. 

The HRHS is constrained by the essentially planar configuration of 
the dual keel space station to move along a direction that is parallel 
to the either the y or a axis of the space station. The position and 
velocity of the HRHS for a given portion of this motion are effectively 
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scalars. The rate and angle of the solar panel displacement are also 
scalars because the rotation of the solar panels is always about the Y 
axis of the space station. 



Motion of MRMS alont? trusses In Y-Z plane only 
Figure 54. Constraints on HRMS motion 


If the system state vector is defined to be a combination of. the 
angular rate of the entire space station, u>, the state of the solar pan- 
els: 11, e, and the state of the HRMS; r^, v^, then a coupled set of 
state equations can be derived by solving the previous three sets of 
equations of motion, eqn's (168) (174) (179) , for the rate of change of 
the state. The equations of motion of the HRHS and solar panels with 
the constraints on their motion added are needed. 

For motion of the HRHS parallel to the vehicle Y-axis, the equations 
of motion for the MRMS become; 
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d^Tj/dt^ » 0 

(F„)y - (a^^dr^/dt) y- ((d«/dt)>‘rT)y - 
0 

And for motion of the NRKS parallel to the vehicle Z-axis, the 

equations of motion become: 

d^rj/dt^ ■ 0 

0 

(Fjj - (Zw^dr^/dt) ((d»/dt)*rT), - 

For the rotation of the solar panels about the Y axis of the vehicle 
alone, the equations of motion become: 

dO/dt - J*’C^ 0 

(Me) - (w^cjn) - (cen^j- jn*)c^«) - («*cjcM - (cjcWdt)„ 
0 (182) 

The notation ()^ and (), refers to the y and z component of the col- 
umn matrix that results from each expression, respectively. By substi- 
tuting these equations for the motion of the solar panels and the MRHS 
into the equations of motion for the entire space station, (168) the 
equations of motion for the core of the vehicle, including the con- 
straints on the NRNS and solar panel motion can be derived. There are 
two sets of these equations, one for motion of the HRHS parallel to the 
Z-ax is and one for motion parallel to the Y-axis. They are, respective- 
ly? 



(180) 
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(Y-motion of HRMS) 


H - (i"d«/dt)^ + + u^cjn)^ + (cCn^J - jn*)cT„ 

(i'Wdt)y + ^ + U^cjn)^ + (C(n*j - jn^)cT« 

(I"du,/dt), + 4 («*cjft), 4 (CCft^J - Jft>‘)cTo 


0 

(«*cjn) , 
0 


(c(n>‘j- jn’‘)cT«)y - («*cjc^w)y - (cjc^d«/dt)y 


4 Kr/ 


4 K U^'r/dr^/dt) J4 K (r/w’^dr^/dt) ^ + K (dr^/dt^w^rJ ^ 
(w^rT^dr^/dt) (rT*<i>’‘dr.^/dt) (dr^/dtVr,) „ 

(«i)’‘r.^*dr.^/dt) j (rT.*«*drT/dt) ^ (drj/dt’'«*r j) ^ 


(Fjy - (2«*drT/dt)y - ((d«/dt) Vy- («’‘«*rT)y 


(2-n»otion of NRNS) 

H - (I-dw/dt)^ 4 U*!"**)^ 4 (»>«CJn)^ 4 (C(ft*J - JO*) CM. 

(I'Wdt) 4 (w*i"«) 4 Mcjft) 4 (c(n*j - jn*)cMv 

(I!'d«/dt), 4 4 (w^cjn), 4 (C(n*j - jn*)cMj 


4 0 

(Mp) - («>‘C 

0 


i/ 

cjn) - (ccn^j- jn»<)cM - (w^cjcM.- (cjc^dw/dt) 


4 Kr/ 


4 K (<tt*‘ry*dr.p/dt) ^ 4 K (rj^w^dr^/dt) ^ 4 K (dr^/dt^w^r^) ^ 

(w’^rT’^dr^/dt) (r/«’‘drT/dt) (dry/dt^r^) . 

(w^ry^dry/dt) j (ry^w^dry/dt) ^ (dry/dt*‘«i)*ry) ^ 


(FJ, - (2«*dry/dt), - ((d»/dt)Vy),- MwVy), 


These equations can now be solved for dw/dt, in terms of a new 


inertia' matrix, I''', and a new 'moment* vector, M' 
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(Iii" - Kr2=) 

+ Kr^rj) 
T •* 


(I " 
VJ.12 


(1,2" + 1,3" 

{I22- - J2 - Kr, 2 ) 

V V 



du)^/ctt 



dwj/dt 

“ 


dcD^/dt 

1 


H/ 

"2' 

"s’ 


(Z motion of HRHS) 


where 1 ,2,3 ■ x,y,z 


’d„" - Kr32) (1,3- + Kr3r,) 


dw^/dt 


r«ii 

1 12" (^22" ■ *^ 2 ^ ^23'' 


dw2/dt 

m 

"2’ 

(1,3" + Kr,r3) I23" (I33" - Kr,2) 


du>2/dt 


"3' 

— — ^ 

<- J 




I' " (Y motion of HRHS) 


The state equations for Y motion of the HRHS are; 


d«/dt • I" 'M' 

dft/dt • J'’Cf[(Hp)y- (»’‘CJft)y- (C(n"j - Jn*)CMy- U^CJCMy 
de/dt ■ n 

dVj/dt - [(Fjy “ (Zw^dr^/dt) y - ( (d»/dt) ’‘r,.) y - («i>’‘<i>*‘rT) y] 

dr^/dt ■ Vj 

The state equations for Z motion of the HRHS are; 


(185) 


(186) 


(CJC^dw/dt) y] 
(187) 
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du)/dt - I" 'M' 

dn/dt » J‘''C^[(Mp)y- («xcjn)y- (C(n’'j - w) (w*CJCVy- (CJC^du/dDy] 

de/dt - n ( 188 ) 

dv^/dt - [(Fjj - (aw^dr^/dOj - ( (du,/dt) 
dr^/dt ■ Vf 

These state equations are coupled and nonlinear, and also have time 
varying coefficients. They have been left in the full nonlinear form as 
a linearization would have put constraints upon the magnitudes of the 
rates and attitudes for which the equations are valid. A linearization 
is usually employed when the goal is a stability analysis of the system. 

In this case we are assuming that whatever control system is employed 
will be able to deal with any passive instability of the space station. 

The equations are in a form which is suitable for numerical integration 
in order to determine the time history of the state; u, ft, e, Vy, r^. 

Once this is known, the momentum buildup on the space station due to the 
modeled disturbance torques may be found. 
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